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I. INTRODUCTION 



One of the main challenges in hadron physics is to understand within Quantum Chromo- 
dynamics (QCD) the structure of the nucleon and its excited states. Through their couplings 
with the meson-baryon continuum, these excited states may be identified with the nucleon 
resonances (A^*) in various meson-baryon reactions. Ideally, one would like to study the 
A^* structure by analyzing the meson-baryon reaction data completely within QCD. This 
however is far from our reach. As a compromise, one can analyze meson-baryon reaction 
data by developing reaction models within which both the internal structure of hadrons 
and the reaction mechanisms are modeled using theoretical guidances deduced from our 
current understanding of QCD. Of course, the separation between the reaction mechanisms 
and the internal structure of hadrons is by no means well defined theoretically. Thus the 
model dependence of the interpretations of the extracted A^* parameters is unavoidable. 
This is not very satisfactory, but seems to be the only option we have at the present time. 
With this understanding, many models have been developed in recent years to analyze the 
very extensive electromagnetic meson production data which have been obtained by intense 
experimental efforts at Jefferson Laboratory (JLab), MIT-Bates, LEGS, Mainz, GRAAL, 
Bonn, and SPring-8. For a comprehensive review of the present experimental and theoretical 
status we refer the reader to Ref. l|. 

In this work we focus on recent measurements of pion electroproduction at energies near 
the AX1232) resonance, and in particular the data at < 2.00 (GeV/c)^ obtained at 

ze these data 
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JLab [2I, |3|, MIT-Bates [4] and Mainz . One possible approach to ana 

is to apply the djTiamical model developed by Sato and Lee (SL) in Refs. 



The SL 



model, as well as the unitary isobar model MAID 10(], and the Dubna-Mainz- Taipei (DMT) 
model Uj have been used extensively to extract the 7A^ A form factors from the data. 



The SL model gives reasonably good descriptions of the very extensive data accumulated 
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since 200 
LEGS 
from J 
data I2 



using parameters obtained by analyzing the pion photoproduction data from 
and Mainz [l^, and the pion electroproduction data at = 2.8,4.0 (GeV/c)^ 
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4,^ 



n^this paper, we would like to re-visit this model and to include the recent 
in a comprehensive analysis to extract more precisely the 7A^ A form 
factors. Furthermore, we explore the theoretical interpretations of the extracted 7A^ A 
form factors in terms of two relativistic quark model calculations \v\ and a recent 
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Lattice QCD calculation [l8l |. 

In section II, we recall the essential ingredients of the SL model and present a refinement 
of the model from improvement of fits to the vrA^ scattering data. Section III is devoted to 
defining explicitly the 'yN A transition form factors which are the focus of this work. The 
extraction of A form factors from electroproduction data is discussed in section IV. 

We then discuss possible theoretical interpretations of the extracted A form factors 

in section V. In section VI, we give a summary and discuss possible future developments. 



II. THE SL MODEL AND ttN SCATTERING 



The essential feature of the SL model is to have a consistent description of both the nN 
scattering and the electromagnetic pion production reactions. This is achieved by applying a 



unitary transformation method 



8l, 



191] to derive an effective Hamiltonian from the interaction 



Lagranrians with A^, A, tt, p ,uj, and photon fields. The details of the model can be seen 
in Ref. |8|, l9|]. For our present purposes, it__is sufficient to just recall the expression of the 
7A^ TiN amplitude TT,N,-yN given in Ref. 



where the non-resonant amplitude is defined by 

tnN,'rN{E) = [tT,N,TvN{E)GT,N{E) + l]vT,N,-yN 

with the non-resonant vriV scattering amplitude defined by 



(2.1) 



(2.2) 



(2.3) 



Here Vj^n^-kN is a vrA^ potential and Gt^n{E) is a vrA^ propagator with relativistic kinetic 
energy operators. 

The resonant amplitude is given by 



r{E) 



(2.4) 



where the dressed vertex interactions are defined, to the first order in electromagnetic cou- 
pling, by 

{E)G^n{E) + 1]T^a,.n, (2.5) 

rA,7Ar = rA,7Af + rA,7rAfG'^Ar(-£^)^7rAf,7Af • (2.6) 



In the above equations, T-^^^a and rA,7Ar are the bare vertex interactions. The energy shift 
Sa(-E') is calculated from 

(2.7) 

Here we remark that the A propagator of Eq. fl2.4p is the mathematical consequence of 
the Hamiltonian formulation of the reactions. It differs from the covariant forms, mainly 
used in tree-diagram models and extensively discussed in the literature on possible off-shell 



effects 



20 



2l| . The investigations on which covariant form of the A propagator is more 
consistent with the general principles of quantum field theory are continuing [22 1. 
By using Eqs. fl2.2p and (12.31) . one can write Eq. (12.61) as 

rA,7Af = rA,7Af + 5TA,'yN , (2.8) 

where 

'^rA,7Af = rA,7rArG'^Ar(-E')l'7rAr,77V (2.9) 

with 

^A,nN = rA,7rAf[l + C^A^ (-£')^7rAf,7rAf] • (2.10) 

We will discuss the dynamical content of the dressed vertices Eqs. (l2.8l) - (I2.10I) in section III. 

Within the considered model the non-resonant ttN potential Vj^n^ttN contains four mech- 
anisms: the direct and crossed nucleon terms, p exchange, and crossed A term. Apart from 
the standard irNN coupling constant, the model is determined by three coupling constants 
QpNNQp-K-K and Hp for the p-exchange, /jtAta for the vrA^ — * A vertex. Each interaction vertex 
is regularized with a dipole form factor F{k) = (A^/(A^ + /c^))^ where k is the pion momen- 
tum associated with the vertex. The model thus has additional three parameters: A^attv for 
the nNN vertex, A,rAfA for the vrA^A vertex and Ap for both pNN and pyrvr vertices of the 
p-exchange term. These six parameters are adjusted along with the bare mass rriA of the A 
to fit the empirical vrA^ scattering phase shifts. We note that this model has considerably 
fewer parameters than the other nN mode 



s used in the study of pion electroproduction cal- 
culations. For example, the nN model [23(] used in the Dubna-Mainz- Taipei (DMT) model 
of pion electroproduction has 16 parameters. Thus we do not attempt to fit the data at 
energies above Tl ~ 250 MeV (invariant mass W ~ 1280 MeV). At higher energies, we 



expect that the couphng with the nnN channel as well as the tails of higher mass A^*, such 
as the Pii(1440) and 5*11(1535) resonances, must be included in a realistic description of the 
ttN scattering data. An extension of the SL model to higher energies has been developed 



recently in Ref. 



The original SL model (model L in Ref. [8]) provides a quite good description of the 
ttN scattering phase-shifts up to Tl ^ 250 MeV except in the P13 partial wave. This can 
be seen from comparing the data with the solid curves in Fig. [1] This discrepancy raised 
some questions concerning the stability of the results with respect to the vrA^ fit. In this 
work, we find that the fit to P13 phase shifts can be considerably improved while retaining 
a similar quality of fit in the remaining partial waves. This is done by weighting the data 
in P13 slightly more heavily in the fit. The resulting fits are the dotted curves in Fig. [H 
The fitted parameters for the new irN model (SL2) and the original SL model in Ref. sl 
are compared in table [H The remaining discrepancy in the 5*31 channel at energies above 

~ 150 MeV perhaps can be resolved if we add additional mechanisms or use a different 
form factor parameterization. We have explored these possibilities and found that the model 
SL2 is the best we can have within the SL model. 

It should be noted here that due to the lack of energy independent solutions at lower 
energies in the P13 and P31 partial waves, we have instead used the SAID |25| energy de- 
pendent solutions and assigned those points with ^ 5% errors in the fits. A more correct 
procedure is to fit the original vrA^ observables. This however is not pursued in this work. 



III. THE 7*iV ^ A FORM FACTORS 

We now discuss the 7A^ — > A vertex interaction defined by Eq. (12.81) . It has a bare 
vertex rA,7Ar and a term 6TA,-yN which is determined by the non-resonant interaction Vt,n,'^n 
and a dressed vrA^ — A vertex f a.-ttAt- In the SL model, v^N^-yN contains the usual direct- 
and crossed-A^ terms, vr-exchange, contact term, crossed-A term, and p- and cj-exchanges. 
Thus the dressed 7A^ A contains the meson loops illustrated in Fig. [2l We note that 
the mechanisms illustrated in Fig. [2] are similar to those in the calculations of the current 
matrix element < A|j^„e^|A^ > within a hadron model where A^ and A contain a pion 
cloud: \B >= Z^^/'^[\Bq > +Ci\Bon > +...]. Such calculations can be found, for example. 



in Ref. [26| using the cloudy bag model. Thus the term 6T^^.yN is called the meson cloud 
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T^^ (MeV) 



FIG. 1: Phase shifts of vrA^ elastic scattering up to Tl =250 MeV. Sohd and dotted stand for model 
SL and SL2 respectively. Data, L2t,2J, are from the energy independent SAID 25| analysis plus 8 



points from their energy dependent solution for the P13 and P31 partial waves at lower energies. 

contribution. Note that this definition is not shared by the DMT modeL 

In a consistent derivation using the unitary transformation method within the SL model, 
all strong interaction parameters except the uNN vertex in v-j^i^^^n have been determined 
in fitting the vrA^ scattering data. The ujNN vertex coupling is determined by using the 
pion photoproduction data at Q"^ = 0, as will be further discussed later. By using the 
previous works on the F^nn{Q'^), F^t,.„{Q^), and -F^7rv(<5^) (with V = p,u!) form factors, the 
non- resonant interaction v^n,ttN is then fixed at any Q^, as explained in Ref. j9j. Since the 
dressed vertex rA,7rAr in Eg. (12.91) has also been fixed in fitting the nN scattering data, the 
meson cloud term ST/^^^n can be predicted at any Q^. With Eg. (12.81) . the bare form factor 
rA,7Af is then the only unknown, and can be extracted from the pion electroproduction data. 
This strategy of the SL model is also used in this work. Ideally, one should fit the vrA^ and 
electroproduction data simultaneously to minimize the errors in extraction. But this will be 
worthwhile only when the polarization data, such as those obtained by Kelly et al. {27], are 
also included in the fits. Our effort in this direction will be reported elsewhere. 

One of the main objectives in this work is to explore the theoretical interpretations of the 
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Model 


Parameter 


SL2 


SL [8] 


9 pN N 9 p-K-K 


37.133 


38.4329 


kp 


2.8496 


1.825 




3.6177 


3.2551 




7.4639 


6.2305 




3.2132 


3.29 


Bare mA 


1296.9 


1299.07 




1.1664 


1.207 


9ujNN 


11.85 


11.82 


Ga/(0) 


1.86 


1.85 


Ge{0) 


0.0251 


0.025 


TABLE I: Parameters of the two considered 


models. The notations 


are explained in the text. 



extracted bare —>■ A vertex. We thus need to define this quantity precisely in the rest 
of this section. The bare vertex Tj\^^n of Eg. (12. 6 p is parameterized in the form developed 



by Jones and Scadron [28|. With the normalization < k\k' >= 6{k — k') for the plane wave 
states and < 4'b\4'b' >= Sb,b' for B = N and bare A states, we can write in the rest frame 
of the A and for the photon momentum q = {u, q) 

x[M^,^,A,A.(g)GM(g') + E^^^,,^,MGEm 
+Cm,,,x,xMGcm] (3.1) 

where and A^v are the helicities of the initial photon and nucleon, mj^ is the z-component 
of the A spin, m^^ and rritj^ denote the isospin components, and 



EN{q)+mN 1 3(mA + m7v) 



(27r)3/2\/ 2Er,{q) v^4m^(^7v(g) + m^) 
with 

Mmj^,x^x^iq) = < nijJiS X q- ex^\\N > (3.3) 




FIG. 2: Graphic representation of the dressed — > A vertex defined by Eq.(2.8)-(2.10). 

-Em,■^,A^A^r(g) = < nijJ^S ■ ex^a ■ q + S ■ qa ■ ex^\XN > (3.4) 
Cm,^,A^AAr('?) = — < rrijJ^S ■ qa ■ qtQ\\N > (3.5) 

where e = A/47r/137, photon polarization vector is defined by e-ti = -^{x ±iy), and eo = 
for = ±1, eo = 1 and eo = for the scalar component = 0. The transition spin S is 
defined by < ii\mis.\Sm\iN'fnN >=< Ja^aUn^^n^ >■ 

The form factors GMiQ"^), Ge{Q^), and GciQ"^) describe magnetic Ml, Electric E2, and 
Coulomb C2 transitions. Choosing the photon direction q in the z-direction, the above 



matrix elements are related to the helicity amplitudes defined by Particle Data Group [29 1 
(PDG) 

M/2{Q'^) = B < rrij^ = 3/2, mt^ = mjj^|rA,7Ar(g)|A^ = +lAAr = -1/2, m^^^ > (3.6) 
^i/2(<5^) = B < nij^ = 1/2, rut^ = mt^|rA,77v(g)|A^ = +1, A^v = 1/2, rut^ > (3.7) 
Si/2{Q^) = B < rrij^ = 1/2, mt^ = mt^\TA,'rN{q)\>^^ = 0,Xn = 1/2, nit^ > (3.8) 

with 
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where 



= ^ (3.10) 

With q chosen in the z-direction, we can easily evaluate Eqs. fl3.3l) - fl3.5p . Eqs. fl3.1l) and (13. 6p - 
(13. 8p then lead to the following explicit relations 



A„2{Q') = -^[GMiQ') + GEiQ')] , (3.11) 
A. 



Ai/2{Qn = -^[GMiQn - ^GEm\ , (3.12) 

Sy2m = -^^GciQ')- (3.13) 
v2mA 



with 



A = ^ '^"^ (3.14) 

2mN Y rriNK^ 1 + Q^/ {mN + 

The helicity amplitudes Eqs. (l3.6p - (I3.8I) are most often calculated in hadron structure 
calculations. Eqs. (13. 110 - 03. 130 then allow comparisons of such calculations with the bare 
form factors G^jiQ'^), Ge{Q'^), GciQ"^) extracted from our analyses. These form factor 
are presumably sensitive to the short-range interquark interactions. They perhaps mainly 
contain information about the quark wavefunctions of the N and A within the constituent 
quark model, as will be discussed in section V. 



In Ref. [9|], it was found that the pion photoproduction data |12l . Il3l | and the electro- 
production data at = 2.8,4.0 GeV^ [l^ can be fitted well by using the following naive 
parameterization 

= ( l + QV0.71(GeV/c)j '^ + '"?^> ^"P*-"'?'' ("^' 

where x = M,E,G, a = 0.154 (GeV/c)^ and b = 0.166 (GeV/c)^. The strengths at = q 
are found to be Gm{0) = 1.85, Ge{0) = 0.025. The value of Gc{0) was fixed using the long 
wavelength limit 

2 

G'c(O) = -4^^^^^^^Gij(0) (3.16) 



mA — m 



N) 



Note that this ansatz assumes an identical dependence for all three 'bare' electro- 
magnetic couplings. There is no theoretical justification for this simple choice. Thus it 
is not clear whether the discrepancies between the predictions from the SL model (using 
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the parameterization Eqs. fl3.15l) - fl3.16p and the data accumulated since 2001 reflect infor- 
mation about the bare couphngs or about deficiencies in the model description of the vrA^ 
rescattering process. With more data in a wide range of Q^, we will abandon in this work 
the parameterization Eqs. fl3.15H3.16l) and directly extract these bare couplings by fitting the 
data at each individual point. 

The dressed form factor rA,7Ar has the same symmetry property of the bare vertex defined 
above. Thus it can be expanded in the same form of Eq. fl3.ip . We denote the dressed 
quantities by Gm[Q'^\ Ge^Q"^), GciQ"^)- The corresponding helicity amplitudes Ax can 
also be calculated by using the same relations Eqs.f l3.111 -( I3.13I) . In Ref. j^, it is shown that 
the dressed ratios can also be calculated from the imaginary parts of Mi+, Ei+ and 5*1+ 
amplitudes of pion electroproduction 

- _ Ge _ ImEi+ 
Gm ImMi+ 
P _ \q\ Gc _ lmSi+ 
2mA Gm ImMi4. 

It is common to define G%.j for the Ml transition form factor which is related to our 
dressed form factor by 

Gm(Q') = M \ „ (3-19) 

where = 115 MeV is used in extracting the data from M^^ amplitude of pion electro- 
production amplitude and T^^ = 93 MeV from the SL model. 



IV. EXTRACTION OF 7iV ^ A FORM FACTORS 

With the refined vrA^ model SL2, we first re-analyze the vr^ photoproduction data. Here 
we need to also tune the less well determined parameters of the a;-exchange ttN 



non-resonant interaction. Following the procedure of Ref. js], we adjust the uNN^ coupling 
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Ij of 



constant guiNN and Ga/IO) and (^^(O) of the bare A^-A form factor to fit the data 
differential cross section and photon asymmetry of the irN reaction. The uNN 

form factor is assumed to be the same as that of pNN which has been fixed in the fit to vrA^ 
data. The quality of the fit is comparable to what was achieved in Ref. [sl, and needs not 
be discussed here. The resulting values of 6*^/(0) and G_b(0) and g^jNN are also listed in the 
last three lines of Table [H 
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0.2 



Q' fGeV" 



FIG. 3: Magnetic dipole transition form factor G\j for 7*A^ A(1232), normalized to the proton 
dipole form factor Gd{Q'^) = + Q^/A^]^ with = 0.71 (GeV/c)^. Experimental points are 
analyses of inclusive data (Q) from pre-1990 experiments at DESY and SLAG [30(1 and recent 



1 



3li. Solid 



exclusive p{e,e'p)7r° data (■) from BATES MAMI 0,3 and JLAB 
curve is from the dressed calculation of this work using the parameterization of Eqs. (|3.15p - (|3.16p . 
The dotted curve is obtained when the meson cloud effect, defined by Eq.(2.9) is turned off. 



Our next step is to extract the A^-A form factors by fitting all of the available pion elec- 
troproduction data at energies close to the A position. As a reference point for our analysis, 
we first use the parameterization Eqs. fl3.15p - fl3.16p of the SL model to make predictions. 
Our results, using that simple parameterization, for the dressed Ml form factors G\.[{Q'^) 
and the ratios Rem and Rsm are the solid curves shown in Figs. [3] and HI The data shown 
here are from the analyses performed by the experimental groups using mainly the unitary 
isobar models 10|, l33|] . The data from our analyses, shown later in FiglHl are not included 



in Figs. [3]and|4l The presence of two different extracted values for the same Q^, in these 
and following figures, is due to the existence of two inconsistent data. 

It is clear that the resulting dressed G\j{Q'^) (solid curve) agree well with the available 
empirical values. In the same figure, we also show the result (dashed curve) which is obtained 
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0.5 



Q' (GeV') 



FIG. 4: The ratios Rem and i?^;^/ (defined by Eqs.(3.17)-(3.18)) calculated from the dynamical 
model using the parameterization Eq.(3.15)-(3.16) (solid line) . The data are from^ecent exper- 



iments at BATES: □ 



MAMI: 



3i, JLAB/CLAS: T [3], ▲ [2] and 



[511, O Ipj, A 13, 

JLAB/Hall A: • [27]. CLAS points from Ref. [3] were obtained from the average of UIM and SL 
model estimates in Table 2. 



by setting the meson cloud effects, STj^^^n defined by Eq. fl2.9p . to zero. Clearly, the meson 
cloud effects are important in the low region and gradually diminish as increases. 
The predicted Rsm (lower part of Fig. H]) also agree well with the data except at the two 
lowest Q2 values Q"^ =0.06, 0.127 (GeV/c)^ from the analyses by the MIT-Bates and Mainz 
groups. 

Here we remark on the data presented in Figs. [SH Because the experiments so far do 
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FIG. 5: The dressed — > A form factors (solid curves) calculated using the parameterization 
EqsJ 



r 



3J5|) 


-(3 


1, A 





not have complete measurements including all polarization observables, the extraction of 
multipole amplitudes from the data needs some constraints imposed by making theoretical 
assumptions. The most common practice is to analyze the data starting from the amplitudes 
generated from the K-matrix isobar model MAID. We note that the data from the MIT- 
Bates and Mainz groups do not have the full coverage of angles and thus their extracted 
values of Rem and Rsm perhaps depend very much on the dynamical content in MAID. 
The results from JLab at > 0.16 (GeV/c)^ are also analyzed using a unitarized isobar 
model (UIM) of Aznauryan 331] . However the JLAB data cover almost the whole angular 
region and hence the fitted results are closer to a full partial wave analysis. 

The ratios Rem shown in the upper part of Fig. H] clearly indicate that the parameter- 
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ization Eg. (13. 151) is not valid for Ge{Q'^)- This can be seen more clearly in Fig. [5] where 
we compare each A^-A form factor with the empirical values. To extract more precisely 
the A^-A form factors, we therefore abandon the parameterization Eqs.f l3.l5]) - fl3.16p and 
perform fits to the available differential cross section data at each by adjusting the 
values of the bare form factors. In Fig. El we show some results (solid curves) from our 
fits at = 0.06,0.127,0.2,0.9,1.45 {GeV/cf. We see that the fits to the differential 
cross section data are better than predictions (dotted curves) using the parameterization 
Eqs. (l3.15|) - (l3.16p . The most visible fit improvement occurs in the longitudinal-transverse 
interference cross section cttl at = 0.06,0.127 (GeV/c)^ 

The resulting bare form factors are shown in Fig. [71 We now turn to discussing how these 
data on the bare A can be interpreted theoretically. 



V. THEORETICAL INTERPRETATIONS 



To explore possible theoretical interpretations of the extracted A parameters, it is useful 
to first consider the electromagnetic pion production reaction within the general reaction 



theory 3J]. Like any reaction involving composite systems, such as the atomic and nuclear 
reactions, its amplitude has a non-resonant part and a resonant part T = t + . Qual- 
itatively speaking, the non-resonant amplitude t is due to the fast process through some 
direct particle exchange mechanisms, and the resonant amplitude t^ is due to the time- 
delayed process that the incoming particles lose their identities and an unstable system is 
formed and then decay into various final states. The unitarity condition ImT = T'^T implies 
that t and t^ are not independent from each other. In the model considered in this work, 
this can be seen in Eqs. (l2.4p for t^ which contain the non- resonant amplitude t through 
Eqs.([23])-(I22D- Thus the extracted dressed form factors G^f (Q^), GdQ'^), Gc{Q^) of the 
resonant amplitude can only be compared with the hadron structure calculations of cur- 
rent matrix element < A|j^^ ■ e^|A^ > which contain the meson loops illustrated in Fig. [2l 
Furthermore, these mesonic effects are dynamically identical to those in the non-resonant 
amplitude t^^jy^^^ = (1 + t7riv,7rAr)f7rAr,7Ar which is equally important in describing the data 
of reaction cross sections. We emphasize here that this close relation between the resonant 
and non-resonant amplitude is not specific to the formulation considered in this work, but 
is the consequence of a very general unitarity condition. 
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p(e,e'p)TT° W=1232 MeV 0'=0.127GeV' 




:0S COS -3,: COS W (GEV) 
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1-1 c ■ - 1 
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FIG. 6: Fits to experimental p{e, e'p)TT^ structure functions. Solid lines are from the fits with 
the bare form factors GuiQ"^), Ge{Q^) and GciQ"^) adjusted at each . The dashed curves are 
from the calculations using the parameterization Eqs.(3.15)-(3.16). Data are from MAMI at 
g2 = 0.06 GeV2, BATES [4] at = o.l27 GeV^, CLAS |^ (W^ = 1220 MeV) and MAMI 
{W = 1221 MeV) at = o.2 GeV^ and CLAS |3] at = o.9, 1.45 GeV^. 
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FIG. 7: Bare form factors for the 7A'^ — > A transition as a function of . The points have been 
obtained by performing individual fits for each value to the correspondingpion electroproduction 



data (BATES: u\^\, MAMI: ■ [4 O 



6fi, A 



% and JLAB/CLAS: T A [0].) The dashed curves 
are from the front form quark model calculations of Refs. [isl. [3^. The dotted curves are from the 
instant form quark model calculations of Ref. [16]. 

Clearly, it is a much more difficult task to understand the extracted dressed G'j\/(Q^), 
Ge{Q'^), Gc{Q^) form factors within QCD. A simpler approach is to assume that the meson 
loops included in our model are the dominant mesonic effects resulting from the spontaneous 
breaking of chiral symmetry. It is then reasonable to compare our bare form factors G'a/(Q^), 
Ge{Q^), GciQ"^) with hadron structure calculations in terms of only the constituent quark 
degrees of freedom. It is instructive to first consider the naive s-wave non-relativistic quark 
model within which the magnetic Ml transitions have a well known relation /xa+p/ /^p = 
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2v/2/3 where 

^//p = (p,m,^ = l\J2^^cTi{z)\p,ms^ = I) (5.1) 
^^f^A+p = (A+,m,^ = |l5^^fTi(2;)|p,m,^ = i) (5.2) 

Here /ip is the proton magnetic moment with physical value /i^^^ = 1 + ~ 2.77, and rriq 
is the constituent quark mass. From the above relation and the definition Eq. fl3.ip . one 
observe that the magnetic Ml form factor of A at = can be directly calculated 

from the proton magnetic moment 



Gm{0) = [V2Gp{0)] 



2{EN{q)+mN) 



2^^^^) - 0.84/ip (5.3) 



ENiq) + ruN 



3(mA + mjv) 

where q = (m^ — m\)/2m/^ ~ 260 MeV/c. We thus observe that our bare value G'm(O) = 
1.85 can be understood in terms of constituent quark degrees of freedom only when we 
assume that the nucleon magnetic moment within the constituent quark model should be 
about ~ 20% less than the empirical value. The difference is from the meson loops, similar to 
that illustrated in Fig. [21 Such a mesonic correction is indeed close to what has been found in 
the cloudy bag model 35| , and is needed to explain the discrepancy between the data and the 
predictions from a covariant model based on the Dyson- Schwinger Equations {sgI. On the 
other hand, our extracted bare E2 transition form factor G'e(O) cannot be understood within 
the non-relativistic constituent quark model. With the tensor force within the conventional 
one-gluon-exchange, the estimated E2 transition of ^ A is known to be negligibly small 
compared with the value calculated from our value Ge{^) = —0.025. 

We now discuss constituent quark model calculations of A form factors. Since the 

considered here is not very small, it is essential to perform calculations relativistically. We 



focus on two recent results from Capstick and Keister |l^ and Julia-Diaz and Riska 16|. 
Both are within the framework of relativistic quantum mechanics outlined by Dirac 37|. 
There are three possible approaches within this framework: instant form, front form, and 
point form. The generators of instant form and front form can be defined within relativistic 
quantum field theory and hence are more closely related to the dynamical reaction model 
considered in this work. The calculation of Ref. ISj is performed within the front form. The 
wavefunctions are expanded up to = 6 harmonic oscillator basis states in a variational 
calculation using a relativitized Hamiltonian ^] which has a one-gluon exchange (OGE) 
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short- distance interaction and Y-shaped string confinement. The tensor interactions and 
spin-orbit interactions are included in both the OGE and confining potentials. All of the 
potentials are smeared with some quark size, and momentum dependence of the potentials 
is parameterized in order to simulate off-shell effects. A quark form factor is included in 
this calculation js^. We assume that this quark form factor is to account for the finite size 
of constituent quarks, not phenomenologically to include the meson loops as illustrated in 



Fig. [21 This conjecture is o: 
The calculation of Ref. 



course debatable and should be explored in future. 



16j starts with wavefunctions parameterized as 0o(-P) = ^(1 + 



|^)~" with P = y 2{k + q"^), where k and q are two intrinsic momenta of the three-quark 
system, and is a normalization factor. The parameters b, a and the constituent quark mass 
mq are adjusted to fit the nucleon electromagnetic form factors. With a phenomenological 
addition of d-state wavefunction, the same parameters are used to predict the 7A^ — > A form 
factors. Here we only consider their calculation within the instant form, mainly because it 
give the best fits to the data. 

The comparisons of the predictions of Refs. SSj and 16[ are given in Fig. [71 Clearly their 
results are only in qualitative agreements with the extracted bare form factors. In particular, 
both calculations fail to reproduce the very rapid drop of GEiQ"^) from the photon point 
Q"^ = to = 0.06. Both calculations predict a sign change of GciQ"^) in disagreement 
with the extracted values. 

In the fits (solid curves) shown in Fig. [6] for extracting the bare form factors, we of 
course also extract the dressed form factors. Their differences, which reflect the meson 
cloud effects deflned by Eq.f l2.9p . can be seen by comparing the solid squares and triangles 
in Fig. [HI Obviously, pion cloud effects on Ge{Q^) and Gc{Q^) are very pronounced at low 
Q^, as predicted in Ref. [^]. In the same flgure, we also display the results from a recent 
Lattice QCD (LQCD) calculations [l8!]. 

We now discuss the results from LQCD calculations. Currently LQCD can only be 
performed reliably with very large quark mass. There are two possible ways to use the 
extracted 7A^ A form factors to test these results. The flrst one is to conjecture that the 
meson-loop contributions in the quenched calculations is suppressed in the calculations with 
large quark mass and hence their results should be compared directly with our bare values. 
The second one is to apply the chiral extrapolation to get results in the physical region with 
correct current quark masses. There are two problems in using the chiral extrapolation. 
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Rem{%) 






Rsm{%) 






UIM 


SL 


SL2 


UIM 


SL 


SL2 


0.16 


-1.94(0.13) 


-2.45(0.2) 


-2.57(0.2) 


-4.64(0.19) 


-4.44(0.35) 


-4.36(0.35) 


0.20 


-1.68(0.18) 


-2.21(0.2) 


-2.31(0.2) 


-4.62(0.18) 


-4.23(0.35) 


-4.14(0.35) 


0.24 


-2.14(0.14) 


-2.70(0.2) 


-2.76(0.2) 


-4.60(0.28) 


-4.32(0.35) 


-4.21(0.35) 


0.28 


-1.69(0.27) 


-1.99(0.2) 


-2.07(0.2) 


-5.50(0.31) 


-5.08(0.35) 


-4.97(0.35) 


0.32 


-1.59(0.17) 


-2.29(0.2) 


-2.35(0.2) 


-5.71(0.33) 


-4.87(0.35) 


-4.75(0.35) 


0.36 


-1.52(0.27) 


-1.80(0.2) 


-1.82(0.2) 


-5.79(0.43) 


-4.76(0.35) 


-4.56(0.35) 



TABLE II: Extracted values of ^2/Afl ratio Rem and C2/M1 ratio Rsm 



5i+/Mi+ at = 
of the p{e, e'p)?:^ 



0.16—0.36 GeV'^ from analysis of preliminary results from a CLAS measurement 
reaction. Methods used are Unitary Isobar Model (UIM) and the SL and SL2 models described in 
this work. Errors are statistical only. 

First, it is only valid for low Q^, although it has been used in a rather high region. 
Second, there are higher order corrections on the commonly used chiral extrapolation, which 
have not been under control. The uncertainties due to this problem have been discussed by 



Pascalutsa and Vanderhaeghen |40| . Thus it is not clear what to conclude from Fig. [8] for 
the results from LQCD of Ref. ij. Further investigations are clearly needed. 

Finally, we present our determined dressed Rem and Rsm in the low region where very 
large meson cloud effects have been identified in Fig. [HI Our results are listed in table HTl and 
compared with the values determined using the unitary isobar model (UIM). The difference 
between our values and that from the UIM reflect some mo del- dependence in the extraction. 



VI. DISCUSSION AND SUMMARY 

Within the dynamical model of Refs. j^, we have performed a comprehensive analysis of 
recent pion electroproduction data to extract more precisely the — > A form factors. It is 
found that the predictions from the original SL model are not changed much by improving its 
fit to the vrA^ phase shifts in P13 channel. The fits to the very extensive vr*^ electroproduction 
data near the A position can be improved when the simple parameterization Eqs. fl3.15p - 
(I3.16P is abandoned and the strengths of the bare form factors GmIQ^), ^^((5^), and GciQ'^) 
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FIG. 8: The extracted A form factors. Dark sq uar es (triangles ) are the dressed (bare) 

values. Open squares are the lattice QCD calculation of Ref. [18\. 

are adjusted at each where the data are available. 

In using our results, shown in Fig. [H] and table [ITl to test hadron structure calculations, 
it is important to note that the differences between the dressed and bare A form 

factors are from the meson loops illustrated in Fig.[2l We emphasize here that the parameters 
associated with these loops are the same as that in the non-resonant amplitude t7rAf,7Af and 
tT^N,TTN which are important components in describing both the vrA^ scattering phase shifts 
and pion electroproduction data. Thus the meson cloud effect identified in this analysis is 
well constrained by the reaction data directly. Such a direct connection with reaction data 
is very difficult to achieve in the current hadron model calculations, in particular the Lattice 
QCD. 

As an exploratory step, we have considered two relativistic constituent quark model 
calculations of —>■ A form factors. The predictions from both models differ significantly 
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from the extracted dressed values, in particular in the low region. As seen in Fig. [TJ their 
results for GMiQ"^) and Ge{Q^) qualitatively reproduce our bare values both in shapes and 
magnitudes. This seems to suggest that the hadron calculations without meson degrees of 
freedom could make contact with the reaction data through the bare from factors extracted 
from a dynamical analysis such as the one given in this work. The origin of the differences in 
the Coulomb form factor GciQ"^) seen in Fig. [7] is not clear. It will be interesting to explore 
whether the parameters within both models can be refined to remove these discrepancies 
and also improve their results for the bare Gm(Q^) and GEiQ"^)- 

It is still difficult to interpret the extracted form factors in terms the current LQCD 
calculations. The discrepancies shown in Fig. [S] indicate the need of a better understanding 
of the chiral extrapolation used in obtaining those results as well as a significant improvement 
in LQCD calculations. 

To end this paper, we mention that the non-resonant amplitude t7rAr,77V; which is crucial 
in evaluating the meson cloud effects on the A form factors, can be more sensitively 

determined in the energy region away from the A position. While the SL model considered 
in this work can account for the data in the region oi W ~ 1100 — 1300 MeV very well 
in most cases [l|, a more extensive analysis of the data away from the A position must be 
performed in the near future when polarization data will also become more extensive. Our 
effort in this direction to further improve our understanding of the A form factors is 

in progress. 
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